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Abstract 

The approximated energy eigenvalues and the corresponding eigenfunctions of the spherical 
Woods-Saxon effective potential in D dimensions are obtained within the new improved quantiza- 
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I. INTRODUCTION 



In recent years, the analytic bound state solution of the hyperradial Schrodinger equation, 
in any arbitrary spatial dimension (D > 2), with several potential models has received 
increasing attention in the literature [1]. There are only a few potentials for which the 
radial Schrodinger equation can be solved explicitly for all n and I quantum numbers. One 
of these exactly solvable potentials is the Woods-Saxon (WS) potential [2] which has been 
firstly solved within the Schrodinger equation for s-states (/ = 0) by Bose [3,4]. However, 
the three-dimensional radial Schrodinger equation for the spherical WS potential cannot be 
solved analytically for / ^ states because of the centrifugal term ~ r" 2 . This potential was 
used as a major part of nuclear shell model to obtain the nuclear energy level spacing and 
properties of electron distributions in atoms, nuclei and atomic clusters [5-7]. It was used to 
describe the interaction of a nucleon (neutron) with a heavy nucleus [8] and also for optical 
potential model in elastic scattering of some ions with heavy target in low range of energies 
[9]. 

So far, numerous attempts have been developed to calculate the bound-state energies 
of the WS potential in the framework of the Schrodinger equation by means of various 
methods; the Nikiforov-Uvarov method [10-12] and the shape invariance and Hamiltonian 
hierarchy [13]. However, the solutions in [10,11] are at most valid for Rq = 0, in which the 
WS potential turns out to become the Rosen-Morse-type potential. The two-component 
approach to the one-dimensional Dirac equation with WS potential is applied to obtain 
the scattering and bound-state solutions [14]. The exact solution of the relativistic Dirac 
equation was obtained for a single particle with spin and pseudospin symmetry moving in a 
central WS potential using two-component spinors [15-17]. An approach has been developed 
to obtain the bound state solutions of the scattering of a Klein-Gordon particle by a WS 
potential [18]. Recently, the approximated eigenenergies and eigenfunctions of the Dirac 
equation for the Woods-Saxon potential and a tensor potential with the arbitrary spin-orbit 
coupling quantum number k under pseudospin and spin symmetry have been obtained [19]. 

Recently, Ma and Xu have proposed an improved quantization rule (IQR) and shown its 
power in calculating the energy levels of all bound states for some solvable quantum systems 
[20,21]. The method has been shown to be effective for calculating the bound state solutions 
of the Schrodinger and Dirac wave equations with a spherically symmetric potential [1,22- 
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26]. So far, it has been applied, with great success, to study a great number of potentials 
like the rotating Morse [22,23], the Kratzer-type [24], the trigonometric Rosen-Morse [25], 
the hyperbolic and the second Poschl- Teller-like potentials [26] and the Hulthen potential 
[1] and so forth. Very recently, Gu and Sun [1] have extended the application of the IQR 
to the solution of the D-dimensional Schrodinger equation with the Hulthen potential for 
I ^ using the usual approximation to deal with the centrifugal term [27-29]. 

In this paper, we present a new systematical approach to solve the Schrodinger equa- 
tion in D-dimensions with WS potential for it's energy eigenvalues and the corresponding 
eigenfunctions by means of the IQR method using the Pekeris approximation scheme [30- 
32] to deal with the centrifugal term. This alternative approach has recently shown its 
accuracy in calculating the analytic spectrum of the Hulthen potential for I ^ [1]. Fur- 
ther, we give analytic tests using energy calculations for interdimensional degeneracy, i.e., 
(n, l,D) — > (n, I ± 1, D =p 2) corresponding to the confined D = 2 — 4 dimensional Woods- 
Saxon potential. 

This paper is organized as follows. In Sec. 2, the IQR method is reviewed and extended 
to any arbitrary dimension (D > 2). In Sec. 3, the D-dimensional (D > 2) Schrodinger 
equation is solved by this method with arbitrary /-states to obtain all the energy states of 
the Woods-Saxon potential using the Pekeris approximation. In Sec. 4, we calculate the 
corresponding hyperradial wave functions of the Woods-Saxon potential. In Sec. 5, the 
interdimensional degeneracy is introduced. Finally, some concluding remarks are given in 
Sec. 6. 

II. IMPROVED QUANTIZATION RULE 

A brief outline to the improved quantization rule is presented with an extension to the 
D-dimensional space (D > 2). The details can be found in Refs. [20,21]. The IQR has 
recently been proposed to solve exactly the one-dimensional (ID) Schrodinger equation: 



where the prime denotes the derivative with respect to the variable x. Here /i is the reduced 
mass of the two interacting particles, k(x) is the momentum and V(x) is a piecewise con- 
tinuous real potential function of x. The Schrodinger equation is equivalent to the Riccati 





(1) 



h 
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equation 

<f>\x) + ^xf + k(x) 2 = 0, (2) 

where 4>{x) = ip'(x)/ip(x) is the logarithmic derivative of wave function ifj(x). Due to the 
Sturm-Liouville theorem, the <fi{x) decreases monotonically with respect to x between two 
turning points, where E > V(x). Specifically, as x increases across a node of the wave 
function ip(x), (f>(x) decreases to — oo, jumps to +00, and then decreases again. 

Moreover, Ma and Xu [20,21] have generalized this exact quantization rule to the three- 
dimensional (3-D) radial Schrodinger equation with spherically symmetric potential by sim- 
ply making the replacements x — > r and V(x) — > V c s(r): 

tb r B 

jk(r)dr = N* + fk> (rMdr, k(r) = ^E^EMR , (3) 

where ta and tb are two turning points determined from the relation E nt i = V c g(r), N = n+1 
is the number of nodes of 0(r) in the region E n j > V c s(r) and it is larger by one than the 
number of nodes of wave function ip(r). The first term Nn is the contribution from the 
nodes of the logarithmic derivative of wave function, and the second term in (3) is called 
the quantum correction. It is found that, for all well-known exactly solvable quantum 
systems, this quantum correction is independent of the number of nodes of wave function 
of the system. This means that it is enough to consider the ground state in calculating the 
quantum correction, i.e., 

TB 

Qc = [ k \r)^P-dr = nu. (4) 

rA 

The quantization rule still holds for Schrodinger equation with spherically symmetric poten- 
tial in D dimensions. In what follows, we shall employ this method to solve the Schrodinger 
equation in D-dimensions with WS potential using the Pekeris approximation to deal with 
the centrifugal term. 



III. EIGENVALUES OF THE WOODS-SAXON POTENTIAL 

The D dimensional Schrodinger equation with spherically symmetric potential V(r) for 
arbitrary /-state takes the form 

+ V(r) - E n A Vv, m (r, Sl D ) = 0, (5) 



where the representation of the Laplacian operator V 2 D , in spherical coordinates, is 

V 2 = — + {D - l) 9 - l ( l + D ~ 2 \ (6) 
D dr 2 r dr r 2 

and 

n D ) = rl> nil (r)YT(n D ), ^(r) = r^ D -^ 2 u{r), (7) 

where Y™ (&,£>) is the hyperspherical harmonics. The wave functions ip n ,i,m(r, Qd) belong to 
the energy eigenvalues E ni i and V(r) is taken, in the present study, as the Woods-Saxon po- 
tential in the configuration space and r represents the Z)-dimensional intermolecular distance 

(s*0 ■ 

Further, substituting Eqs. (6) and (7) into Eq. (5) yield the wave equation satisfying the 
radial wave function u(r) in a simple analogy to the 2D and 3D radial Schrodinger equation 

d 2 u(r) 2u ri _, . ., . . „ . . 

+ f [E n ,i - V eff (r)} u(r) = 0, (8) 

where E n j is the bound state energy of the system and V e ff(r) is the deformed Woods-Saxon 
[2] effective potential in D dimensions defined by 

_ (r-flp) 

1 + ge « ' 

with the parameter 

A = ( A2 -^ h \ A = 21 + D - 2. (10) 
8/x 

where Rq = tqA 1 ^ is the nuclear radius with tq = 1.25 fm (may vary by as much as 0.2 
fm depending on the specific nuclide) and A the atomic mass number of protons (Z) plus 
number of neutrons (N) in target nucleus, Vq ~ 50 MeV (having, dimension of energy) is 
the potential well depth and a ps 0.5 — 0.6 fm is the length representing the surface thickness 
that usually adjusted to the experimental values of ionization energies [12]. In the case of 
negative eigenenergies (i.e., when E n j e [—50 MeV, 0]) we have the well-known bound-states 
problem while in the case of positive eigenenergies (i.e., when E n j G [0, 1000 MeV}) we have 
the well-known resonance problem. It should be noted that the deformation parameter q is 
a real parameter and can be taken equal to 1 and —1 for the Woods-Saxon and the Hulthen 
potentials, respectively in the calculations. The radial wave function u(r) satisfying Eq. (8) 
should be normalizable and finite near r = and r — > oo for the bound-state solutions. 
The wave equation (8) with the Woods-Saxon potential is an exactly solvable problem for 



/ = (s-wave) [3,4,15], however, it cannot be solved analytically for A 7^ because of 
the orbital coupling term Ar™ 2 . Therefore, to solve Eq. (8) analytically, we must use the 
Pekeris approximation [30-32] to deal with the centrifugal term. By inserting the conversions 
x = (r — Rq) /Rq and a = Ro/a, where r G (0, 00) 4i£ (—1, 00), the following exponential 
form can be used instead of the centrifugal term 



1 



Rl 



-2ax 



d + di 



+ d 2 



(11) 



which can be expanded around the minimum point r ps Rq (or x = 0) only up to the second 



order as 



Rl 



, di d 2 a . , , . a 2 , 2 
d + — + — - - (di + d 2 ) a: + y^" 2X 



(12) 



where d , di and d 2 are coupling constant parameters. It is worth to note that the above 
expansion is valid for low rotational energy states and q — 1. The expression r~ 2 can be also 
expanded around x = up to the second order term as 

1 1 



Rl 



(l-2x + 3x 2 ) 



(13) 



Comparing Eqs. (12) and (13), we obtain the constant parameters for the Woods-Saxon 
potential as follows 

, , 4 12 , 8 48 , 48 

d = 1 h — , di = 5-, d 2 = — , (14) 



a* 



a' 



However, these constant parameters for the Hulthen potential take values as follows 



do = dl = 1' ^ 2 = 1- 



(15) 



Further, by defining 



Rq 



1+ 



D-3 



(16) 



where the scaling parameter Q is introduced with the purpose of obtaining solution for the 
Hulthen potential as well. So we take Q — 1 for the Woods-Saxon potential but Q = a and 
Rq = 1 for the Hulthen potential. 
Equation (8) can be rewritten as 



d 2 u(x) 2fiRl 



dx 2 



E n ,i + Vq 



1 + qe~ c 



-2ax 



- b ( d + di 



l + qe- 



+ d 2 



'1 + qe~ ax Y 



u[x) 



= 0, 
(17) 
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where n and I signify the radial and orbital angular quantum numbers, respectively. 

We now study this quantum system through the improved exact quantization rule. At 
first, we introduce a new variable 

V = lT^' Tx = ~ ay(1 ~ qV) > Tr = -R- y(1 ~ qV) ' (18) 
where y G (0, x + qea )■ Overmore, the turning points i/a and db are determined by solving 
V eff(y) = bd 2 y 2 + (bdi - V ) y + bd = E n>l as follows: 

VA=^f-^- 7^r^(V -bd 1 ) 2 + 4bd 2 (E nl -bd ), (19a) 



2bd 2 2d 2 2bd 2 

y*=^ 2 -^^dj { y«- bd ^ 2 + m ^ - hd ^ ( 19b ) 

with the properties 



VA+VB=^-±VAVB= ~ (Era - bd ) . (20) 
The momentum k(y) between two turning points is expressed as 



= ^ L \[bd 2 \/{yB -y)(y-y A ), (21) 

dk (y) _ via ^ ( fvKEv _ .fyEK) , (22) 



dy 2h V V y~VA V vb -y 

The Riccati equation (2) now becomes 

- ^y(l - qy)^^ = ~% [E - bd 2 y 2 + (V - bd 1 ) y - bd ] - My) 2 , (23) 
R dy h 

having the only possible first order polynomial solution satisfying 

Mr) = ciy + c 2 , ^1 = ~c lV (l - qy), Cl > 0. (24) 
ar Kq 

where we have used 0o( r ) = 0o(y)- Substituting </>o(y) into Eq. (23), one has the ground 
state wave function and energy eigenvalue solutions 

MV) = ^rny + - jg^ ifr) , 

m = -2 (l ± Jl + f^^l^'l , m > 0, (25) 



8^ Q 2 8 2 d 2 
W a 2 q 2 



for the Woods- Saxon 



The new quantum number is chosen as m = — § ^1 — ^Jl 
potential and m = — | ^1 + \Jl + ^pp^ j for the Hulthen potential. After a lengthy 
algebra but straightforward, we can calculate the integral of the quantum correction (4) 
based on the ground state as 



/ y/2fj,d 2 5 Q m 

Vc = -7r — 1 

\ n. g a q 



1 



(26) 



The integral of the momentum k(r) in the quantization rule (3) is calculated as 

vb i *b 



fj z ( r ) dr= yfifi /^2^o r I y/(y-VA)(yB -y) | y% - yA ) (vb - y) 

J h \1 q 2 a J \ y ' 1 



qy 



dy 



= 7T- 



2 A« #c 



{\[^ + ~~ En ' 1 ~~ ^ (f° 



// a \ \j <•/- ' v V Q ¥ J Q 

Using the relations (26) and (27), the improved quantization rule (3) turn out to be 



(27) 



7T- 



H a 



[\[^ + ~ En ' 1 ~ ^ 
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Vn 



7r n 



^2jid 2 5 Q m 



(28) 



h q a q 

Thus, based on the improved quantization rule we can finally obtain the approximated energy 
levels E n i valid for both the Woods-Saxon and the Hulthen potentials in the D-dimensions, 



n,l 



— "0 



2[x a 2 



x 



« 2 (2» + 1 =F V 1 + lw^) f W + d.) - 



4i? 2 



Q 2 ^ 2 J 

„2„2 «2 



(29) 



i?^2n + l Tv /l , 

where n, I — 0, 1, 2, • • • . Further, by setting Q = 1 and the deformation parameter q = 1 in 
the above equation, we get the energy spectrum of the Woods-Saxon potential as 



E 



(D) 



5 2 J ftV 

rdn X 



R, 



2 a 



2fi 



2n + l- 



8 M s 2 j 



"I 2 



2m il 



2 ^T/ 



4a 2 



(2n + l-0Tf|S) 



(30) 



In three-dimensions (D = 3), it can be reduced to the form 



2 ^2 



(D) _ h 2 l (f + 1) 



2^ 



2// 



(2n + 1 - ^1 + 4/ (/ + !)$) (/ + 1) (d, + d 2 ) - f^ 



"I 2 



4a 2 



(2n + 1 - yfl + 41(1 + 1) J) 

(31) 



which is essentially same as that obtained by other Nikiforov-Uvarov method [33]. 

On the other hand, the case when we set the parameters values do = 1/12, d 2 = d s = 1, 
q — — 1, Q — a, Rq — 1 and V = aZe 2 in Eq. (29), then we obtain the energy spectrum of 
the Hulthen potential: 

2 ' 



F (D) J 1 

n ' 1 2/x I 12 



1 + 



D - 1 



/ + 



D-3 



fiZe 



h 2 (n + l + 2=±) a 



(32) 



which is similar to the one obtained in Ref. [1] (if d = in the usual approximation). 



IV. EIGENFUNCTIONS 

We are now in the position to study the corresponding eigenfunction of this quantum 
system for completeness. The Riccati equation of the relation (8) is 

0V) = -^[^-V; // (r)]-0(r) 2 , (33) 

where 

0(r) = (34) 

Based on 

u(r) = ef r ^ dr = e~ ar W^)* {z)dz , (35) 
and using Eq. (24), we can easily calculate the eigenfunction of the ground state as 

u (r) = N (e ~j (l + qe ~ J , e > 0, v > 1, (36) 

where 

e n=0 = (bdo-Eo), v = ^1 - ^1 + ^^Zj , (37) 

and iVo is the normalization constant. 

Let us find the eigenfunction for any quantum state n other than ground state. By 
considering the boundary conditions 

„ when r — » 00, 

y={ , (38) 

1 when r — > 0, 
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with u(y) — > 0, we may define a more general radial eigenfunctions, valid for any quantum 
number n, of the form: 

u(y) = y^(l + qyYf(y), y = e"^, e n> , > 0, z/ > 1. (39) 

Substituting Eq. (36) into Eq. (8) leads to the following hypergeometric equation 

y (1 + qy) f"(y) + [1 + 2e n ,, + q(l + 2e n , x + 2v) y] f'(y) 

^a 2 (6di - Vb) - ^ (1 + 2e n ,,) /(y) = 0, (40) 
with the following conditions: 



ei tl + ^a 2 (E n j - bd ) = 0, (41a) 
q 2 v(y - 1) - %a 2 bd 2 = 0, (41b) 



which leads to the following results 



, D-l\f 1 D-3\Q 2 d 2/i 0rn . . 

^ = «\/ ( I + -g-J ^ + -j- J - ^o 2 ^, (42a) 



,= i(l ± 0,C=^l + 4(z + ^l)(z + ^)^4, (42b) 

where + and — are taken for the Hulthen and the Woods-Saxon potentials, respectively 
The solution of the Eq. (40) is then 

F{y) - 2Fl (A, B, C, y) - T{A)T{B) ^ f^Tk) W (43) 

where 

A = e n j - qv - VU = -n, n = 0, 1, 2, • • • , 

B = e n>l -qu + VU, (44) 
C = 1 + 2e nh 

where 2 iq (A, B;C;y) is the hypergeometric function. Now, we may write down the radial 
wave functions (34) as 

u n ,i(r) = Nni [e » J (^1 + ge « J 2 Fi ^-ra,n + 2e n) i + g£ - g; 1 + 2e n>J ;e « J 

(45) 
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If now, we set n = in Eq. (45), then we can easily obtain Eq. (36). Finally, the 
unnormalized total wave functions are obtained as 

x 2 F, (-n,n + 2e nil + q(-q;l + 2e n x,e ~j Y™^). (46) 

which is identical to Eq. (42) of Ref. [33] when D — 3. Thus, the Jacobi polynomials can 
be expressed in terms of the hypergeometric functions [34] 

Pl A > B \1 - 2s) = r( "** + ,f 2F 1 (-n, n + A + B + 1; A + 1; x). (47) 
n!l (1 + A) 

It is worth noting that the hypergeometric function 2 F\ {A, B;C;x) is a special case of the 
generalized hypergeometric function [34] 



p F q (a u a 2 , • ■ ■ , a p , fi u fr, ■ • • , fi q] x) = £ ^H^*'^ ^ ^ 
where the Pochhammer symbol is defined by (y) k = T(y + k)/Y{y). 



V. INTERDIMENSIONAL DEGENERACY 

From Eq. (30), it can be seen that two interdimensional states are degenerate whenever 
[35] 

(n, /, D) -+ (n, I ± 1, D =f 2) =► i^I?, (49) 

i.e, 

"I 2 



5 2 , r fl 



2^2 



n,l -*n,l±l ~ R 2"° 2/i 



(2n + 1 - v/TTfJ^) fjf § ( dl + d 2 ) - %V 



(50) 

Thus, a knowledge of E^ for D = 2 and £> = 3 provides the information necessary to find 
for other higher dimensions. 
For example, = E^\ = E^\ = E^l. This is the same transformational invariance 
described for bound states of free atoms and molecules [36-38] and demonstrates the existence 
of interdimensional degeneracies among states of the confined Hulthen potential. 
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VI. CONCLUSIONS 



By using the improved quantization rule, we have given an alternative method to obtain 
the approximated energy eigenvalues and eigenfunctions of the Schrodinger equation in D- 
dimensions with the Woods-Saxon potential for all /-states within the Pekeris approximation 
to deal with centrifugal term. We emphasize that the expressions obtained for the energy 
eigenvalues and eigenfunctions are valid for all real values of Rq. However, the solutions 
provided in Refs. [10,11] are at most valid for R = (Rosen- Morse- type potential and not 
the Woods-Saxon potential). The advantage of this method is that it gives the eigenvalues 
through the calculation of two integrations (26) and (27) and solving the resulting algebraic 
equation. First, we can easily obtain the quantum correction by only considering the solution 
of the ground state of quantum system since it is independent of the number of nodes of 
wave function for exactly solvable quantum system. Second, the wave functions have also 
been obtained by solving the Riccati equation. The general expressions obtained for the 
energy eigenvalues and wave functions can be easily reduced to the Woodes-Saxon and the 
Hulthen potentials, the three-dimensional space (D — 3), s-wave (I = 0) cases. The method 
presented here is a systematic one, simple, practical and powerful than the other known 
methods. Finally, the simplicity of the method motivates us it to solve the Dirac equation 
for this WS potential to see the resulting relativistic effects. 
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Appendix A: Integral Formulas 

The following integral formulas are useful during the calculation of the momentum integral 
and the quantum correction term: 




(Al) 
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1 



zdr 



IT 



r A 



r V( r — r A)ijB — r) yfrxHs' 

dr = 7T. 



y/(r- r A )(r B - r) 



TB 

J ^ \/Jr- taHtb - r)dr = n 



-jir A + r B ) - \/r A r B 



(A2) 



(A3) 



(A4) 



rB 

f- 

J (a 



dr 



7T 



TA 



(a + br) y/(r - r A )(r B - r) yj{a + br A )(a + br B ) 



--, r B >r A > 0. (A5) 
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